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Abstract 

We employ density-functional-theory calculations to analyze the interactions between oxygen 
interstitial atoms and (a)-type screw dislocations ((a) = AbA)) in a titanium, based on inves¬ 
tigations of generalized stacking fault (GSF) energies, dislocation core structures and the strain 
field of an oxygen interstitial. Compared with the substitutional atoms such as Al, there is a large 
repulsive interaction between oxygen interstitials at the octahedral site and an (a)-type screw dis¬ 
location core, since the volume for the interstitial atom between Ti lattice sites is largely reduced 
in the core. Differential displacement maps and distributions of interstitial volume surrounding the 
dislocation core show that the interaction with an oxygen interstitial is very short ranged (~ 5 A) 
and directional, mainly located on the prismatic plane where the majority of lattice displacements 
in the (a) screw dislocation core occurs. At long distances the interactions mediated by the strain 
from the interstitial oxygen and the dislocation vanish. We show that the large repulsive interac¬ 
tion at short range can push an oxygen atom from its original interstitial site to a new interstitial 
site on the basal plane in the dislocation core. The strong repulsive interaction and the mechanical 
shuffle of oxygen interstitials can be expected to have a large impact on dislocation dynamics, and 
can be used to explain the strong strengthening effects of oxygen impurities in a titanium observed 
experimentally. 
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I. INTRODUCTION 


Solid solution alloying and impurity atoms change the strength, plasticity and toughness 
of metallic alloys. Such effects are extremely strong for oxygen and nitrogen impurities 


in Ti and its alloys|l|, which have many important structural applications, especially in 
aviation, aerospace and marine industries due to their high specific strength, excellent high- 


temperature mechanical performance and corrosion resistance^. For example, even an 
increment of 0.1 wt% oxygen impurities in commercial purity a-Ti can significantly increase 
its yield strength dramatically, but with a concomitant reduction in toughness [2], Recently 
in-situ transmission electron microscopy (TEM) compression tests were performed on Ti 
nanopillar samples with different oxygen concentrations [ll . The results show that only 0.3 
wt% O can increase the yield strength to 2.5 GPa for ce-Ti . Further characterization 
and analyses indicate that such high strength mainly results from the interactions between 
oxygen interstitial atoms and screw dislocations with Burgers vector (a) = that 

glide on the prismatic planes {llOOjjsJ. These recent findings have motivated the present 
computational work, aimed at further understanding the nature of the interactions between 
oxygen interstitial atoms and (a)-type screw dislocation in a-Ti. 

Solid solution strengthening effects are generally regarded to originate from two types 
of mechanisms, commonly referred to as elastic interactions and chemical interactions[1, 4], 
The former is due to the long-ranged interactions between the strain fields of solute atoms 


and those due to dislocation lines 


, 6]. The latter can change the dislocation core structures 


and their mobility by varying the restoring force due to lattice slip; this type of interaction is 


rort-ranged and can be related to the so-called generalized stacking fault (GSF) energy[4, 
)]. The elastic interactions between solid solutions and the screw dislocation lines are 
relatively weak; according to isotropic linear elasticity theory, screw dislocations generate 


purely shear strain fields, which wou 


d lead to zero elastic interaction with solute atoms 


that generate only volumetric strains 10]- Another type of elastic interactions is due to the 


change of elastic modulus because of solute atoms[l]. However, such effect is not significant 
for 0.3 wt% O in a-Ti lattice jll]. Thus, the high strengthening effects for oxygen in a-Ti 
more likely result from the chemical effects. 

Using first-principles calculations, chemical effects underlying solid-solution strengthening 
can be investigated both indirectly and directly. Indirectly, first-principles calculations can 
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be used to investigate the change of the GSF energy due to the presence of solute atoms 
on the slip plane, and the resulting effects on dislocation core! structures can be estimated 
based on semidiscrete variational Peierls—Nabarro model jl2. 13]. Directly, the interaction 
between solute atoms and the core structure of a straight dislocation can be calculated either 
by the combination of lattice displacements distribution inside the core structure of pure 
metal and the change of GSF due to the solute atom( 4 ], or by inserting the solute into the 


core structure in first-principles calculations of the core 


14 -In |. 


Short-ranged solute-dislocation chemical interactions can be classified into two categories 
depending on the lattice site for the solute atom. For substitutional solutes, replacing a 
matrix atom on a regular lattice site, the effects one the dislocation core may be relatively 
small as long as the solute has similar atomic radius compared with the matrix atom. Thus, 
it is accurate to estimate the chemical interaction based on lattice displacements inside the 
core structure of pure metals 15, 16]. On the other, the interaction between a dislocation 
core and an interstitial atom may be quite large especially when the latter have relatively 
large atomic radii; in such a case the interaction energy will depend on the interstitial 
volume, which could be quite different in the dislocation core than the bulk lattice, due 
to the large lattice displacements and distortions. For example, previously a model was 
proposed based on TEM characterization that (a) screw dislocation is slightly spread on 
the prismatic plane and would also dissociate into partial dislocations with (c)-components 
on Burgers vector Ji3, ljj]; these partial dislocations have to combine together for further 
slip on the prismatic plane (so-called sessile-glissile transformations ); the oxygen impurities 
could impede such re-combination processes to strengthen Ti 18,LL9]. Recently several first- 
principles studies on the a-Ti (a) screw dislocation core structure indeed show multiple 


complex core structures and there are still debates on the most stable conhguration 


2 ( 3 - 22 ], 


To understand the origins of the significant strengthening effects of oxygen interstitials 
in a-Ti, we employ in this work density functional theory (DFT) calculations of oxygen 
impurity effects on the GSF energy along the (1120){1100} slip system, as well as on the 
(a) screw dislocation core structure in a-Ti. We also calculate the lattice distortions due to 
oxygen interstitials in a bulk Ti lattice to estimate interactions with the strain held of screw 
dislocations, considering lattice anisotropy. Both GSF and dislocation core simulations 
show that, compared with the relatively weak strengthening effects of the substitutional 
solute atom Al, the extremely strong effect of oxygen and other large interstitial atoms 
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result from the strong repulsion resulting from the large changes of interstitial volume inside 
the dislocation core. This strong interaction is found to be very directional (mainly on one 
prismatic plane) and short ranged (<~5 A). In addition, the large repulsive interaction from 
the screw dislocation core can mechanically shift the oxygen from its original interstitial site 
to a new interstitial site on the basal plane in the dislocation core. We also discuss how the 
findings of the present atomistic calculations may affect the behavior of dislocation dynamics 
at larger scales. 


II. METHODS 


A. Calculations of Generalized Stacking Fault Energy 

DFT calculations were performed by using the Vienna ab initio simulation package 
(VASP) 23, 124]- All the calculations were performed non-spin-polarized, using the pro¬ 


jector augmented wave (PAW) method 


M 


26], and the Perdew-Burke-Ernzerhof (PBE) 


exchange-correlation functional 27]. There are 4 and 6 valence electrons for Ti (3d 2 4s 2 ) and 
O (2s 2 2p 4 ), respectively. Partial occupancies of eigenstates were determined by a first-order 


28]. The cutoff energy 


Methfessel-Paxton smearing scheme, with a broadening of a = 0.2 eV 
for the plane wave basis was 400 eV. 

As shown in Fig. |Tj (a), the supercells for the GSF energy calculations were constructed 
with a (2 (a) xl (c)) periodicity within the slip plane, with the axaxis along the (1120) Burg¬ 
ers vector and the y -axis along (0001); along the perpendicular z-axis ((1100)) there were 12 
layers of prismatic planes plus 13 A of vacuum. For this superccll, (5x7x1) Monkhorst-Pack 
k-point grids for the Brillouin-zone integration were applied j^J. An interstitial atom can 
be placed at the octahedral site on the interface, as indicated by the dashed line in Fig. 
H] (a), between two prismatic planes where the slip is set to occur. In the hep structure 
there are two atoms in one primitive cell, one atom on type A and the other on the type B 
basal planes; these sites are indicated by blue and red colors, respectively, in what follows. 
There are two different geometrically-distinct types of prismatic interfaces. However, the 
one indicated by the dashed line in Fig. Q] (a) is found to lead to much lower GSF energies, 
so that the other is not considered further in this work. 

To understand solid-solution effects, a substitutional atom can replace a Ti atom on the 
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prismatic plane nearest to the slip interface. Here substitutional atoms were used as a refer¬ 
ence in order to compare the different strengthening effects from interstitial solutes. During 
GSF calculations, a fault vector x was applied by shifting the top half of the supercell rigidly 
along the Burgers vector (x-axis) direction. Under these fault vectors, DFT calculations were 
performed by relaxing all metallic atoms only along the 2 -axis directionand the interstitial 
atom along all three axes. To check the effects of solid solution concentration, supercells 
with different sizes ((2 (a) x2 (c)) and (3 (a) x2 (c))) were employed in calculations of the 
change of the GSF energy when one oxygen interstitial atom is placed on the slip plane (per 
area of the periodic supercell). The GSF energy results with oxygen interstitials depend on 
the initial position of the interstitial atom before DFT relaxations, so different initial oxygen 
positions were investigated. Further details are given in Section IIII1 


B. Calculations of Screw Dislocation Cores 


We performed a series of DFT calculations of Ti screw dislocation cores without/with 
solute atoms. As shown in Fig. [1] (b), the supercell has 12 (c) periodicity along the first 
vector [0001], 6 (a + c ) periodicity along the second [1101] vector, and 1 (a) periodicity along 
[1120] direction, so there are totally 288 Ti atoms in the supercell for pure a-Ti. There 
are two straight screw dislocation lines along the [1120] direction with opposite Burgers 
vectors: 1UM anc [ t 1120 ! ; respectively. Such a dislocation dipole configuration is called a 
quadrupolar stacking, where opposite-sign dislocations are stacked upon one another in 
a periodic dislocation dipole array to decrease the effects of long-range strain fields on 
dislocation core structures 


30 


31]. The dislocation dipo 


e configuration within the superccll 


32]. Therefore, the lattice vector of 


also results in a distortion of the supercell vectors 
the supercell along [1101] is tilted ^ along the Burgers vector direction to reflect a dipole 
composed of compact dislocations separated by one half the dimension of the supercell along 
the (c) direction. The initial displacement fields surrounding the dislocation core were solved 
algebraically using Daw’s approach |3l(, where the inputs are elastic constants of cc-Ti and 
the geometric centers of the dislocation cores. 

In setting up the initial displacement fields, Ti elastic constants were taken from exper¬ 
imental values, and we used two different geometric centers: one is at the octahedral site 
between two basal planes, the other is the center between two Ti atoms on the same basal 
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planes. These initial structures result in slightly different dislocation core structures as ob¬ 


served in previous studies 


20 


22], which will be discussed in the following sections. After 


relaxation to obtain the dislocation core structures in pure a Ti, A1 as a substitutional solute 
and O as an interstitial solute were placed at various positions relative to the dislocation 
core centers. All atoms in the supercell were relaxed freely, but the superccll shape and 
volume were constrained to be fixed. To save computational time, a “soft” PAW potential 
for oxygen was used, so that the cutoff energy for the plane wave basis could be reduced 
to 282.8 eV in all dislocation core calculations. A 1 x 1 x 11 Monkhorst-Pack fc-point grid 
was used for the Brillouin-zone integrations. The convergence criteria for the force was 0.01 
eV/A. All the other details associated with the DFT calculations were the same as described 
for the GSF calculations above. 

C. Calculations of The Strain Field Induced by an Oxygen Interstitial 


The above calculations are used to calculate the interactions between solute atoms and 
screw dislocation cores at short distances ( < 10 A). At longer distances, the interactions 
between point defects and dislocations can be described by elastic theories based on the 

n n n 

stress/strain fields of the dislocation and those of the individual point defect [5|, [6|, [33]. 


The strain field generated by (a) screw dislocation in a-Ti can be obtained analytically by 
using the elastic constants of a-TijlO, 34], The strength of individual point defect can be 


described by 
calculations 


he stress/strain of this point defect generated in a finite superccll from DFT 
6 ]. Thus, we inserted one O interstitial atom at an octahedral site in a perfect 
Ti hep structure, using an orthorhombic superccll (3 (a) along [1120], 2 (c) along [0001], 2 
a/3 a along [1101], respectively), and calculated the stress/strain due to this O impurity by 
using the same DFT methods for the GSF calculations. The O-induced stress was calculated 
by fixing the superccll and relaxing all the inner coordinates, and the stress-free strain was 
obtained by relaxing both the inner coordinates and the superccll shape/size. 
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III. RESULTS AND DISCUSSIONS 


A. Effects of Solute Atoms on Generalized Stacking Fault Energies 


As determined by the lattice structure, the full GSF energy curve for ct-Ti on the prismatic 
plane for slip along (a) is symmetric with respect to the slip distance x = 0.5 a, which 
corresponds to the maximum misfit across the prismatic interface. Thus we only show GSF 
results from x = 0 to x = 0.5 a in Fig. [2] (a). For pure Ti, the GSF energy at x=0.5 a takes 
the value 7 t;(x = 0.5a) = 206.7 mJ/m 2 , and it corresponds to a shallow local minimum, 


consistent with previous studies 


21 ], 


As mentioned in III A1 the GSF energy calculated with an interstitial atom on the slip 
plane was found to depend on the initial position of interstitial solute before relaxation. Since 
the stable site for the interstitial atom is determined by the available interstitial volume, 
we plot in Fig. [2] (b) the distribution of interstitial volume in the a-Ti structure near the 
prismatic interface as the slip distance gradually increases. Here the interstitial volume of 
any point in the 3D lattice is defined as ^ d 3 , where d is the distance between this point 
to its nearest Ti atom. This definition is expected to provide a reasonable estimate of the 
space available to an interstitial atom, since the energy of the solute atom is expected to be 
strongly affected by its shortest Ti-interstitial bond length. We can see that in the perfect 
lattice (x = 0) there are periodic regions with large interstitial volume located in the middle 
between a type-A and a type-B Ti atom on two nearby basal planes. Each of these regions 
corresponds to an octahedral interstitial site in Fig. Q] (a). The 3D atomic structure of 
oxygen at such an interstitial site is plotted in Fig. [2] (c). 

As the slip distance increases, the interstitial volumes for the octahedral sites on the 
slip interface, labeled with a dashed line, gradually decrease. When the slip distance is 
equal to 0.5a, corresponding to the largest misfit across the slip interface, the interstitial 
volume near the original octahedral site on the slip interface is much less than half of the 
value in the perfect lattice. On the other hand, there are regions on the basal plane near 
the solid interface where the interstitial volume increases with increasing slip displacement, 
and they have almost the same interstitial volume as the original octahedral site in the 
perfect structure, when £=0.5 a. As shown in Fig. [2] (d), this large interstitial volume 
corresponds to a new octahedral site on the basal plane when the slip distance close to 0.5a. 
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The same conclusion regarding the transition of interstitial site was described recently in an 


independent study 


35]. 


Based on the above results, we set two paths for the GSF calculations when there is one 
oxygen interstitial atom on the prismatic interface. In Path I the oxygen is located at the 
original octahedral site between two basal planes before DFT relaxations, and oxygen in 
a nearby local minimum in the energy after relaxation. In Path II the oxygen interstitial 
is at the site that corresponds to the octahedral position on the basal plane for large slip 
distances, as illustrated in Fig. |2](d). However, since this second site is not a local minimum 
in the energy for the undistorted Ti structure, the GSF energy for Path II was calculated 
starting from x= 0.5 a with descending slip displacement. 

Figure [2] (a) shows the calculated effect on the GSF energy due to solute atoms. Results 
are plotted for O in both paths I and II described above, as well as for an A1 substitutional 
atom on the prismatic interface. We first focus on the GSF energy with an interstitial O in 
Path I, which should be a natural choice for O starting from zero slip in the reference perfect 
lattice. Here both A1 substitutional and O interstitial atoms increase the GSF energy along 
path I, as compared with results for pure Ti. These results indicate that both solute species 
can strengthen Ti. In addition, the increase in GSF energy due to O along Path I are much 
higher than those from substitutional Al, for all ranges of slip distances. This suggests that 
the solid-solution strengthening effects of Al are weaker than for O. This strong effect of 
oxygen is easy to understand based on the large change of interstitial volume in Fig. [2] (b), 
which indicates that a very large strain is applied to the oxygen interstitial atom as the slip 
displacement increases. On the other hand, for the substitutional atom the change of atomic 
volume and the corresponding strain on the atoms are generally much smaller jj, Q. 

The calculated GSF energy values strongly depend on the solid-solution concentration 
on the interface, as determined by the superccll-size. Thus, it is necessary to obtain a 
quantitative value that is independent of the solid-solution concentration to describe the 
interaction between the dislocation core and an individual solute atom. We thus convert 
the value of the GSF energy, 7 per area, to a value per atom A E = 7 x i s , where A s is the 
interfacial area in the superccll; this quantity is also labeled in the right y-axis in Fig. [2] (a). 
Further, we define the value Eg^(x) as follows: 


e sa( x ) = (7sa(z) - 7Ti(x)) x A s = AEsa(x) - A E Ti (x) 


( 1 ) 
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where 7Ti(^) and 7 sa(^) is the GSF energy for pure Ti and Ti with a solute atom (SA) on 
the interface under the slip distance x, respectively. For supercells with large A s , E'^(x) is 
the energy difference to transfer a solid-solution atom from its equilibrium site in a bulk Ti 
structure to a site on the slip interface with slip distance x. Since in the center of a screw 
dislocation core there should be regions where x = 0.5a, corresponds to the largest lattice 
misfit across the slip plane, E^(x = 0.5a) can be regarded as a reasonable estimate of the 
interaction energy between a single solute atom and the screw dislocation core. 

The results in Fig. [2] (a) show that E^ PatilI (x = 0.5a) = 1.76 eV E^(x = 0.5a) = 
0.23 eV, further confirming the strong repulsive interaction between the oxygen interstitial at 
the original octahedral site and the dislocation core, as compared with the smaller repulsive 
interaction for the substitutional atom. Increasing the supercell size for GSF calculations 
to (2 (a) x2 (c)) and (3 (a) x2 (c)), £'o-Pathi( x = 0.5a) increased to 1.94 eV, so this 
strong repulsive interaction is not a consequence of high solid solution concentration. The 
magnitude of this repulsive interaction for interstitial oxygen is thus much stronger than the 


solid-solu 
0.1 eV 


ion effects associated with substitutional atoms, which are usually on the scale of 


15 


16]. 


The GSF energy associated with an oxygen interstitial atom along Path II decreases as 
the slip increases. In the 2(a) x 1(c) supercell, Path II becomes more favorable than Path I 
when x > 0.25a. This indicates that the energy can be reduced at high slip displacements 
by transferring oxygen from the original octahedral site, between two basal planes in the 
undistorted Ti structure, to the new octahedral site on the nearby basal plane. However, 
such an oxygen shuffle may require overcoming a high energy barrier, as oxygen diffusion in 
the bulk Ti structure is known to have barriers as high as 2 eV 36]. Thus, we estimate the 
barrier for the oxygen shuffle for two paths at fixed slip distance. Because the shuffle distance 
for oxygen is small and the shuffle path is simple, the barrier for oxygen shuffling at constant 
slip displacement can be estimated by linear interpolation of the atomic positions along the 
transition path, where the reaction coordination is the oxygen position along the [0001] axis. 
In these energy-barrier calculations, for each intermediate step the oxygen coordination is 
fixed along the [0001] direction but relaxed along the two other directions, and all metallic 
atoms are only relaxed along the z-axis to fix the slip distance 

The results are shown in Fig. [2] (e). When the slip displacement is 0.50 a, oxygen at 
the octahedral site in Path II near the type B basal plane is 1.62 eV more stable than the 
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site in Path I, and the energy barrier for shuffling between the two sites is 0.42 eV. When 
the slip displacement is 0.25 a, the stability of oxygen in the two paths is almost equal and 
the barrier for oxygen shuffling between the two sites is 0.95 eV. In addition, since there are 
two symmetry-equivalent directions for Path II, for oxygen near type A and type B basal 
planes, respectively, we investigate another oxygen shuffle path from Path I to Path II near 
a type A basal plane when the slip displacement is 0.30 a. The results show that oxygen at 
the octahedral site in Path II near the type A basal plane is 0.49 eV more stable than the 
site in Path I, and the energy barrier for shuffling between the two sites is 0.74 eV. So it can 
be estimated in the real dislocation core, the barrier for oxygen to shuffle between the sites 
in the two different GSF paths should be on the scale of 0.5~1.0 eV, which is comparable 


with the large oxygen diffusion barrier in the bulk Ti structure 


36]. The large scale of the 


interaction energies and the energy barriers for oxygen shuffling (compared with x 300 
K = 0.0259 eV) suggest that oxygen (and most likely other similar large interstitials such 
as nitrogen) are strong dislocation pinning points. The large diffusion barriers also imply 
that the effects of dislocation slip and plastic deformation due to these pinning obstacles are 
likely to be sensitive to temperature and strain rate. 


B. Interactions between Solid Solution Atoms and Screw dislocations 

The above GSF calculations can be used to estimate indirectly the effects of solute atoms 
on dislocation core structures. However, in the real dislocation core, the slip distance and the 
corresponding lattice misfit is not a constant value but has an inhomogeneous distribution. 
In addition, there can be relative displacements between atoms other than those exactly 
located on the slip interfaces. To examine this possibility, we performed direct first-principles 
simulations of screw dislocation cores in a-Ti, and study its interactions with solute atoms, 
as described in the next section. 


1. Dislocation cores in pure a-Ti 


The relaxed configurations of the (a) screw dislocation cores in pure ct-Ti are plotted in 
Figure [3] based on Vitek’s differential displacement (DD) map 37]. In these DD maps, the 
direction and magnitude of each vector between two nearest atoms stand for the relative 
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displacement along the Burgers vector direction ([1120], perpendicular to the plane of the 
paper) in the dislocation conhguration compared with its counterpart in perfect lattice. In 
addition, the integration of the differential displacements surrounding a closed circuit on this 
map gives the total Burgers vector in the region surround this circuit. Here Fig. [31(a) is the 
DD map for the dislocation core conhguration relaxed from the elastic solution where the 
geometric core center is located at the octahedral site between two basal planes; Fig. [3] (b) 
is the map for the dislocation core relaxed from the elastic solution where the geometric core 
center is located on the basal plane. In both maps most of the displacement fields are spread 
on one prismatic plane along the [0001] direction, which is consistent with experimental 
reports that (a) screw dislocation mainly glides on a 2D prismatic plane 3,ll9|. However, in 
Fig. El (a) there are also relatively large displacements on another prismatic plane nearest 


to the glide plane. So Fig. [3] (a) an 


and symmetric core, respectively 


20 


b) are commonly referred to as the asymmetric core 
21], Thus, a slight difference of initial conhguration 


determined by the geometric core center in the elastic solution results in different dislocation 
core structures. In our calculations under the given supercell structure, the asymmetric 
core is about 0.025 eV per Burgers vector more stable than the symmetric core, roughly the 
equivalent of room temperature. Although the energy differences may depend on the external 
strain held that varies with supercell structures, this small energy difference suggests that 
there are multiple possible metastable states of dislocation core configurations, and their 
energy differences are very small so that it is not difficult for the dislocation core structure 
to transform from one to the other 22]. 

There are two atoms in each primitive cell in the hep structure that can be labeled as 
atoms on type A and type B basal planes, denoted by blue and red colors, respectively. We 
plot another type of DD map to calculate the differential displacement vector only between 
the same types of atoms (either type A or type B), which can be used to describe the relative 
displacements between different primitive cells. Compared with DD maps between every 
atom, DD maps based on only type A/B atom can be easily used to identify the location of 
dislocation core by integration of displacement vector surrounding each primitive cell. For 
both the asymmetric core in Fig. [3] (c) and symmetric core in Fig. [3] (d), the integrated 
differential displacement surrounding a primitive cell formed by four nearest type-A basal 
atoms always show that the screw dislocation core can be located in just one primitive cell: 
the integrated displacement of the lattice rectangle surrounding the primitive cell is exactly 
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one Burgers vector (a), but for all other primitive cells the integrated displacements are 
exactly zero. Same behaviors are also found for DD maps based on only type-B basal atoms 
for both cores. It means that even though the dislocation core structure is widely spread, it 
is still not dissociated into two partial dislocations; otherwise there should be at least two 
primitive cells that include finite Burgers vectors. We also increased the size of the supercell 
along the (c) direction by changing the periodicity of the supercell to 18 (c) x 4 (a + c) and 
still did not observe dislocation dissociation. 

Following the GSF calculations in Fig. [2] (b), we plot the distribution of interstitial 
volume surrounding the asymmetric and symmetric core in Fig. [3] (e) and (f), respectively. 
In both figures, the green dashed rectangle is the primitive cell that contains integrated 
displacements equal to a full Burgers vector, as shown in Fig. [3] (c) and (d). Consistent 
with the GSF results, in both dislocation core centers, there is more than a 50% reduction of 
the interstitial volume at the original octahedral site in the dislocation core, which is labeled 
by the dashed circle inside the primitive cell with full Burgers vector in Fig. [3] (e) and (f); 
this is regarded as the dislocation core center in the pure Ti lattice in this paper. In addition, 
there are also large increases of interstitial volume on the two basal planes near this original 
octahedral site, where the new interstitial volumes on the basal planes are almost equal to 
those octahedral sites in the perfect lattice, consistent with the new interstitial site on the 
basal plane when slip = 0.5 a in Fig. [2] (b) and (d). Thus, both the displacement maps 
and the changes of interstitial volumes confirm that there is a finite region in the center of 
the dislocation core with relative displacements very close to 0.5a. According to the GSF 
calculations, the repulsive energy between the dislocation core and an oxygen interstitial 
atom can be much larger than 1.0 eV, and such interaction can cause the oxygen to move 
to the new octahedral sites on the nearby basal planes. 

2. Solute atoms in the dislocation core 

The next step is to study the interactions between solute atoms and (a) screw dislocation 
cores. First, for substitutional atoms, there are two Ti lattice sites on the two basal planes 
closest to the dislocation core center in pure Ti, indicated by dashed circle in Fig. [3] (e) and 
(f). We replace a Ti atom in either of such two sites for both symmetric and asymmetric 
cores in pure Ti with A1 atom and perform DFT relaxation calculations. Two types of 
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DD maps with A1 substitutional atom at such sites on type-B basal planes for the original 
asymmetric core are shown in Figure |4] (a) and (b), where the filled green circle represents 
A1 and the open golden circle is located at the same place of the dislocation core center 
in pure Ti indicated by dashed circle in Fig. |3] (e). There are changes for each individual 
differential displacement vector by comparing Fig. [I] (a) and Fig. [3] (a), The primitive cell 
which contains integrated displacements equal to a full Burgers vector is still located at the 
same position compared with the original DD map for pure Ti (Figure |4] (b) vs. Fig. [3] 
(c)), and there are also no other primitive cells with non-zero integrated displacements in 
Figure d] (b) as found also in pure Ti. This implies that the A1 substitutional atom does 
not change the position of the dislocation core. We also obtain the same conclusion for A1 
at the other basal plane location (type-A basal) and in the symmetric core. The results are 
consistent with the finding that the repulsive interactions between substitutional solutes and 
the dislocation core structures are relatively small (on the scale of 0.1 eV for many alloying 
cases) HQ lbj ] . It should be emphasized that the A1 solute in these calculations has a high 
concentration because of small superccll dimension along (a); nevertheless, it is still found 
to be stable inside the dislocation core. This is consistent with the relatively small increase 
in the calculated GSF energies induced by A1 substitutional atoms relative to pure Ti. 

We also insert an oxygen interstitial into the dislocation core, which is located exactly 
in the position of dashed circle in Fig. |3] (e) and (f). Two types of DD maps with O 
interstitials at such sites in the symmetric core are shown in Figured] (c) and (d), where 
the golden filled circle represents the position of the O atom, inside the golden empty circle 
located at the same place of dashed circle in Fig. |3](f). As shown in Fig0(c), the DD map 
between each atom changes significantly from its original planar shape to a distribution on 
two nearby prismatic planes. Correspondingly, in Figdl (d) we can see that the primitive 
cell with one integrated Burgers vector is pushed away to anther prismatic plane. We find 
similar results for the case of the asymmetric core. These results imply that the repulsive 
energy between oxygen atom and the dislocation core is so large that the dislocation core is 
effectively pushed away by the interstitial solute. This is consistent with the large value of 
-^o-Pathi( x = 0.5a), and the significant reduction of interstitial volume as slip increases in 
the GSF calculations. On the other hand, the DFT relaxation does not induce the oxygen 
solute to move into adjacent sites, such as the new octahedral site on the basal plane in 
the dislocation core; this result is again consistent with the large activation barrier for an 
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oxygen to shuffle to nearby site, as discussed in the previous section. 

Previous DFT simulations with different boundary conditions suggest that Ti [1120] screw 
dislocation core may dissociate into two partial dislocations because of the local minimum in 
GSF located at slip = 0.5 a in Fig. [21(a) 21], However, if the dissociation exists, the region 
between two dissociated dislocation cores should be a stacking fault layer with relative slip = 
0.5 a. According to our GSF calculations, this means that the region with small interstitial 
volume would further grow, and would increase the repulsive interaction between oxygen 
interstitial atoms and the dislocation core further. 


3. Solute atoms near dislocation cores 

To further clarify the magnitude and character of the strong repulsive interactions be¬ 
tween an oxygen interstitial and a dislocation core, we also insert oxygen solutes at different 
interstitial sites. Two types of interstitial positions were investigated. One type is the oc¬ 
tahedral site on the same glide prismatic plane of the dislocation core with distance ±-y, 
±(c) and ±yy, respectively. The other type is on one of the two nearest prismatic planes 
closest to the glide prismatic plane of the dislocation core, and the distance between the core 
octahedral site to the investigated site along [0001] is 0, ±yr, ±(c) and ±^y-, respectively. 
In total, there are 20 octahedral sites investigated for asymmetric and symmetric cores. To 
calculate the interaction energies between oxygen interstitials at those sites and in the dis¬ 
location core, a reference state is also calculated by inserting one oxygen in the supercell at 
a distance >16 A away from the dislocation core, for both asymmetric and symmetric con¬ 
figurations. Based on this reference, the interaction energies between the oxygen interstitial 
and (a) screw dislocation core are calculated. 

The histogram of oxygen solute-dislocation interaction energies for all 40 configurations 
(asymmetric and symmetric) is plotted in Fig. [5] (a). It clearly shows that almost all the 
interaction energies are relatively small, and on the scale of ± 0.06 eV, which is comparable 
to x T = 0.0259 eV, at T = 300 K. So there are relatively weak interactions between 
the oxygen and the (a) screw dislocation core, as long as the oxygen is not located inside 
the region with large lattice displacements. There is one exception, with repulsive energy of 
~ 0.1 eV, corresponding to the case that an oxygen interstitial is y- away from the center 
octahedral site; the corresponding DD map is shown in Fig. ED (b). We can see that the 
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position of the dislocation core center, the primitive cell with one integrated Burgers vector, 
is shifted away from the original position labeled with the golden open circle for pure Ti in 
Fig. [31(f). According to Fig. [31(e) and (f), the interstitial volume of the octahedral sites ^ 
away from the center octahedral site is also largely reduced and their shapes are also largely 
distorted, indicating that there could be also significant repulsive interaction between the 
dislocation core and the oxygen interstitial. Thus, the dislocation core center can also be 
pushed away from the original position as the DD map shown in Fig. [5] (b) indicates. 

Fig. [31(c) and Fig. [31(d) show DD maps between type A atoms for two other typical cases 
in all 40 configurations. The former corresponds to the case that the oxygen interstitial is 
(c) away from the original octahedral center at the same prismatic glide plane; the latter is 
the case where O is at the prismatic plane nearest to the glide plane with 0 distance along 
[0001] direction to the original octahedral center. Apparently, although the introduction of 
the oxygen interstitial changes the structures of dislocation cores to a certain extent, the 
position of the dislocation core center, i.e., the primitive cell with one integrated Burgers 
vector, does not change. This is consistent with the weak interaction energies shown in 
Fig. [31(a) and also the distribution of interstitial volume in Fig. [31(e) and (f). The largest 
changes of interstitial volumes occurs for sites on the same prismatic plane of the dislocation 
core with the range of ~ (c) away from the original octahedral center. Thus, it can be 
concluded that the interaction between oxygen and an (a)-type screw dislocation core is 
very short ranged (~ 5 A) and is very directional, being most pronounced for sites on the 


same prismatic plane as the major glide plane of the (a) screw dislocation il Il9]. 


4- Long-range interactions between oxygen interstitial and (a) screw dislocation 


The weak interaction between an oxygen interstitial and the screw dislocation for sites 
outside the small core region is also consistent with the linear elasticity theory. In isotropic 
elasticity theory, a perfect screw dislocation generates a pure shear stress held, which does 
not interact with a solute atom that creates simp, volumetric stramQ . The interaction 

even when anisotropic elasticity of dislocation 


remains relatively weak (on the scale of 0.1 eV 

fr¬ 
aud the solute atom are taken into account [4j-6, [15j]. To check the effect of anisotropy in the 

present case, we first calculate the stress/strain induced by the occupation of oxygen at an 

octahedral site in the perfect Ti lattice in a octahedral supercell defined in Sec. Ill Cl The 
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results show that, besides the volumetric changes, oxygen at octahedral sites only generate 
shear stress/strain on the (1120) plane perpendicular to [1120]. This result can be easily 
understood by the projections of atomic positions of Ti and O atoms shown in Fig. |T| 
(a). It shows that projections of Ti atoms on (0001) plane (the left part of Fig. [T] (a)) 
have mirror symmetry relative to (1120) and (1100) across the position of the O atom, and 
similar mirror symmetry can also be found on (1100) plane. So the chemical bonds between 
this O atom and all nearby Ti atoms also have such mirror symmetry on these two planes, 
which produce no shear deformation. But on the (1120) plane (the right part of Fig. |T] 
(a)), no mirror symmetry exists relative to (0001) or (1100) across the projection of the O 
position. On the other hand, in a hep structure a straight screw dislocation along [1120] 
generates only shear strain components on (0001) and (1100) without other volumetric or 
shear components [lO, 34]. Thus, by a first-order approximation, there should be no elastic 
interactions between this screw dislocation and O interstitials located at octahedral sites in 
a hep structure. 


5. Mechanical shuffle of oxygen interstitial 

All the above calculations characterize the interaction strength and range of the interac¬ 
tions between an (a)-type screw dislocation and an oxygen interstitial solute located at its 
normal octahedral site, which is the most stable site in the undistorted Ti structure. How¬ 
ever, as shown in Sec. IIII A1 there could be other stable sites in the dislocation core when 
there is large misfit across the slip plane, and these sites cannot be reached automatically in 
the DFT relaxations, due to the large barrier for oxygen to shuffle inside Ti. Thus, instead of 
the original octahedral sites between two basal planes, the oxygen is also inserted at the new 
octahedral site on either of two basal planes in the relaxed configuration of the dislocation 
core, where there is large interstitial volume according to Fig. |3](e) and (f). Both the asym¬ 
metric and symmetric cores are investigated and the relaxed DD maps for the symmetric 
core case are shown in Fig. [6] Similarly to those for A1 substitution in Fig. [4] (b), although 
the internal atomic structures change, the position of the dislocation core center determined 
by the primitive cell with one Burgers vector integrated displacement does not change (this 
is also found for the asymmetric core, not shown here). In addition, compared with oxygen 
far away from the dislocation core, it only costs 0.05~0.10 eV extra energy to insert oxygen 
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into the basal plane of the dislocation core. From the GSF energy calculations in Sec. IIIIA1 
-^o-Path n( x = 0.5a) = 0.14 eV in (2(a) x 1(c)) superccll and converge to 0.10 eV in (3(a) x 
2(c)) supercell. These results suggest the repulsive energies estimated from direct dislocation 
core-interstitial interactions are on the same scale of its approximate value E 1 Q t (x) derived 
from GSF energy calculations. The differences come from the inhomogeneous distribution 
of lattice displacements inside the dislocation core and the lateral interactions between oxy¬ 
gen in the image supercells along (a) direction, which is also close to 0.1 eV according to 


previous DFT calculations[38j. All these structural and energetic results suggest that it is 


likely that oxygen interstitials can exist in a stable form at such new octahedral sites on the 
basal plane in the dislocation core, which is not possible in the undistorted structure of a-Ti 
due to the close-packed structure of basal plane. In fact, such oxygen interstitial structures 
were observed by high-resolution scanning-transmission electron microscopy (HR-STEM) in 
a-TiQ], which is a strong evidence to support the validity of our first-principles simulations. 


IV. SUMMARY AND DISCUSSION 

In summary, we have employed first-principles DFT calculations to investigate the in¬ 
teraction between A1 and O solute atoms and [1120] ( a )-type screw dislocations in a-Ti. 
Both GSF energy and dislocation core calculations show that, when compared with a sub¬ 
stitutional solute atom like Al, a large interstitial atom like O has a much stronger but 
short-ranged and directional repulsive interaction with the (a)-type screw dislocation. This 
strong repulsive interaction for the interstitial case can even lead to the mechanical shifting 
of an oxygen solute from an original octahedral site to a new position on the basal plane in 
the dislocation core. This strong interaction qualitatively explains the large solid solution 
strengthening effects of oxygen to pure Ti found in experiments. While the above calcula¬ 
tions are performed for oxygen interstitials, many of the results correlate with the changes 
in interstitial volumes induced by the dislocation, such that we expect qualitatively similar 
results for large interstitial solutes like C and N, more generally. 

In our simulations, because of the high oxygen concentration of interstitial solutes arising 
from the small superccll size along the (a) direction, the dislocation core can be pushed 
away as a whole straight line when it is close to oxygen at the octahedral site. However, in a 
realistic case when the oxygen concentration is small, other plastic deformation mechanisms 
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are much more energetically favorable. One possible scenario is that only local parts of 
the screw dislocation line near the oxygen interstitial cross slips to the nearby prismatic 
plane. In this case, the segments of dislocation lines that connect the two parts of the screw 
dislocations on two nearby prismatic plane can be quite immobile jj| , which may even act 
as two pinning points as a Frank-Read sourcesjioj]. Another possible mechanism is that 
the dislocation core pushes the oxygen interstitial away from the original octahedral site 
to the new site on the basal plane. Then as the dislocation core sweeps away, the oxygen 
interstitial atom has two possible choices. It can either jump back to their original octahedral 
site if the perfect lattice is almost recovered, or it can be pushed away to next octahedral 
site in front of the dislocation line. Under large plastic deformation where there are many 
screw dislocations gliding on the prismatic plane, this multiple-step oxygen shuffle could be 
a possible mechanism for oxygen interstitial atoms to accumulate in front of certain screw 
dislocation lines, which would further impede the mobility of the dislocations. Thus, large- 
scale simulations, based on the basic interaction mechanisms obtained from this paper, 
should be applied to quantitatively understand the strengthening effects and toughness 
degradations for oxygen and other large interstitial atoms in a-Ti. 
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(a) (b) 

Figure 1: Supercells for (a) GSF and (b) dislocation core calculations. Ti atoms in type A/B basal 
planes are labeled as blue/red color. In (a) the dashed line is the location of slip interface, and 
oxygen atom is labeled as golden color on the interface. In (b) the clockwise and anti-clockwise 
circle stands for screw dislocation core with Burgers vector = -plilAll. 
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Figure 2: (a) GSF for pure Ti and Ti with solid solution atoms on the interfaces. The left y- 
axis is labeled in unit of interfacial energy; the right y-axis is to show the corresponding energy 
difference in the (2 (a) xl (c)) super cell, (b): Distribution of interstitial volume in a-Ti lattice 
near the prismatic interface as the slip distance gradually increases. At any point, the maximum 
value of interstitial volume along [1120] is plotted, (c): Lattice structures of oxygen (golden) at 
octahedral site between type A (blue) and type B (red) basal planes in perfect Ti lattice, (d): 
Lattice structures of oxygen at octahedral site on type B basal plane when slip = 0.5 (a), (e): The 
energy paths for oxygen shuffle from (c) (reaction coordinate = 0) to (d) (reaction coordinate = 
1) at different fixed slip distance. 
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(e) 


(f) 


Figure 3: (a) and (b): DD maps between each Ti atom for asymmetric core (a) and symmetric (b) 
core in pure Ti. The lattice orientation is the same as Fig. |T](b). (c) and (d): DD map between 
only type A (blue) Ti atoms for the same asymmetric core and symmetric core, respectively. Here 
* is located in the primitive cell with the integrated displacement equal to one Burgers vector, 
(e) and (f): Distribution of interstitial volume in the same asymmetric core and symmetric core, 
respectively. Dashed rectangles are located in the same primitive cell with * sign in (c) and (d). 
Dashed circles indicate the region that corresponds to an octahedral site in perfect lattice but with 
large reduction of interstitial volume in the dislocation cores. 
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Figure 4: (a) DD map between each Ti/Al atom when A1 atom (green) is inserted in the asymmetric 
core. The golden circle is located in the same position as the dashed circle in Fig. [3] (e). (b): DD 
map between only type A Ti atoms for (a), (c) DD map between each Ti atom when O atom 
(golden dot located in the golden circle) is inserted in the symmetric core. The circle is located in 
the same position as the dashed circle in Fig. [3] (f). (d): DD map between only type A Ti atoms 
for (c). In both (b) and (d) * is located in the primitive cell with the integrated displacement equal 
to one Burgers vector. 
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Figure 5: (a) Histogram of interaction energies between asymmetric/symmetric dislocation core 
and oxygen atom at various positions, (b)-(d): DD maps between only type A Ti atoms for oxygen 
at several typical positions. The meaning of each symbol is the same as Fig. [I] (d). The golden 
circle is located in the same position as the dashed circle in Fig. [3] (f), which can be defined as 
original dislocation core center, (b): oxygen is ^ away from the original center, (c): oxygen is 
(c) away from the original center, (d): oxygen is at the nearby prismatic plane with zero distance 
along (c) from the original center. 
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(a) (b) 

Figure 6: (a) and (b): DD map between only type A Ti atoms for oxygen at type B (a) and type 
A (b) basal plane near the original dislocation core center in the symmetric core, respectively. The 
meaning of each symbol is the same as Fig. [4] (d). 
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